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Observation of flat band for terahertz spoof plasmon 
in metaUic kagome lattice 
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We study the dispersion relation of a metamaterial composed of metalhc disks and bars arranged 
to have kagome symmetry and find that a plasmonic flat band is formed by the topological nature 
of the kagome lattice. To confirm the flat-band formation, we fabricate the metamaterial and make 
transmission measurements in the terahertz regime. Two bands formed by transmission minima that 
depend on the polarization of the incident terahertz beams are observed. One of the bands corre- 
sponds to the flat band, as confirmed by the fact that the resonant frequency is almost independent 
of the incident angle. 
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Kagome lattices have attracted considerable interest 
from the aspect of geometric frustration in condensed- 
matter physics.^ There is extensive degeneracy of nondis- 
persing resonant modes in a resonator system with 
kagome symmetry. These eigenmodes form a flat band, 
where the resonant frequency of the band is the same 
for all wave vectors in the first Brillouin zone. The flat 
band originates purely from the topology of the lattice 
structure, and it remains flat even when the couplings 
between the adjacent sites are significantly large. Fur- 
thermore, such flat bands can lead to ferromagnetism 
of itinerant fermions,^^ supersolidity for bosons, ^'^ crys- C 
talline ordering,^ and other effects. 

Although the flat-band formation is first expected 
in quantum systems, it can occur in electromagnetic 
systems. The presence of electromagnetic flat bands 
in kagome lattices has already been predicted the- 
oretically in some electromagnetic systems, such as 
two-dimensional photonic crystals^ and metallophotonic 
waveguide networks. ^^ In the flat band, the group ve- 
locity is slowed down in all directions and the effective 
mass of the photons becomes very heavy. It is impor- 
tant to study the flat band in the electromagnetic system 
with kagome symmetry in terms not only of fundamental 
physics, but also from an application standpoint, such 
as slow light; ^^ however, there has been no experimental 
demonstration for the electromagnetic flat band. 

Here, we focus on the flat band for a terahertz (THz) 
plasmonic mode. Although there is no surface plas- 
mon of metals in the THz region, the modes having 
the dispersion relation similar to surface plasmons are 
formed in structured metals, and called spoof surface 
plasmons. ^^^^ In this paper, we theoretically and exper- 
imentally obtain the dispersion relation for a spoof plas- 
mon in the metallic kagome lattice and demonstrate the 
electromagnetic flat band in the THz regime. Numerical 
simulations are also performed to provide confirmation 
of the experiments. 



We introduce kagome-type bar-disk resonators (KB- 
DRs), shown in Fig. 1(a). Metallic disks and narrow bars 
are connected to form a kagome lattice. KBDRs are arti- 
ficially engineered metallic structures, and considered as 
a kind of metamaterial. ^^^^ In KBDRs, electric charge 
stored on each disk temporally oscillates between posi- 
tive and negative values. We discuss the formation of a 
flat band in KBDRs by using a coupled oscillator model. 
We denote the electric potential at the zth disk as (j)i. 
Introducing <I>^ := J (pidt^ we obtain a Lagrangian C as 
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with capacitance C of the disk, inductance L of the bar, 
coefficient of electric induction Cm between nearest disks, 
and adjacency matrix Aij of the kagome lattice, whose 
element is 1 if the ith and jth disks are directly connected 
by a bar for i y^ j; otherwise 0.^^ The first, second, and 
third terms of Eq. (1) represent the electric energy stored 
on disks, mutual electric energy stored between disks, 
and magnetic energy stored around bars, respectively. 
Here, we consider only the nearest mutual couplings. 
Using the Euler-Lagrange equation {d/dt){djC/d^i) — 
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FIG. 1. (Color online) (a) Schematic view of kagome-type 
bar-disk resonators (KBDRs). (b) The fabricated KBDRs on 
a stainless-steel plate with / = 800, d = 10, and r = 145, and 
thickness h = 30 /xm. 
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FIG. 2. (Color online) (a) Dispersion relation of KBDRs for 
/€ = 0. (b) Localized eigenmode of KBDRs. 



dC/d^i = 0, we obtain coupled charge equations as 

qi + cJq (% - ^ Aijqj) + /^ ^ Aijqj = 0, (2) 



with stored charge qi = C ^j a t the ith disk, resonant 
angular frequency cjq = 1/a/ZC, and coupling coefficient 
hz = Cm/C In the frequency domain, we rewrite Eq. (2) 



as 






(3) 



where tildes represent complex amplitudes. Owing to the 
lattice symmetry, we can reduce Eq. (3) to an eigenvalue 
problem for a 3 x 3 matrix and obtain the dispersion 
relation consisting of three bands as 




where F = cos/cy -ai +cos/c|| •a2 + cos/c|| • (ai — 02) with 
wavevector k\\ in the xy plane and unit-lattice vectors 
{ai, 02} shown in Fig. 1(a). A calculated dispersion 
relation is shown in Fig. 2(a) for /^ = 0. It is clear that the 
highest band uj/ujq = ^^6/(1 — 2/^) is flat or independent 
of /c||. It can be seen that the lowest band shows conical 
dispersion near the F point and that the bending middle 
band touches the flat band at the F point. 

The flat band is caused by the interference of spoof 
plasmon propagating in the kagome lattice. The adja- 
cency matrix Aij of the kagome lattice has eigenmodes 
localized at hexagonal sites with an eigenvalue of —2. 
One of them is shown in Fig. 2(b). The number of the 
eigenmodes is equal to the number of hexagons in the 
kagome lattice. The flat band is formed from these de- 
generated localized modes as they are not coupled with 
each other. 



III. EXPERIMENTAL SETUP 



We fabricate KBDRs on a stainless-steel plate. The di- 
mensions depicted in Fig. 1(a) are as follows: period be- 
tween bars / = 800 /im, bar width d = 10 /im, disk radius 
r = 145 /im, and metal thickness h = 30 /im. The size of 
the area patterned KBDRs is 1.1 cmx 1.1 cm. A photomi- 
crograph of a fabricated sample is shown in Fig. 1(b). 

To investigate the dispersion relation experimentally, 
we perform THz time-domain spectroscopy (THz-TDS), 
shown in Fig. 3. A THz emitter and detector (EKSPLA 
Ltd.) with dipole antennas attached to Si lenses are 
used. These antennas are integrated on low-temperature- 
grown GaAs photoconductors, and driven by a femtosec- 
ond fiber laser (F-lOO, IMRA America, Inc.) with a wave- 
length of 810 nm and pulse duration of 120 fs. The THz 
beam is collimated with the Si lens near the emitter. 
The beam radius is about 3.7 mm, which covers a large 
number of KBDRs. The THz electric field E{t) is coher- 
ently measured by the detector. We obtain the trans- 
mission spectrum T{uj) in the frequency domain from 
T{uj) = \Es{uj) / Erei{uj)\'^ , whcrc Es{uj) and £^ref(^) are 
Fourier transformed electric fields with and without the 
sample, respectively. 

In order to obtain the band structure between the F 
point and the M point, the sample is rotated by 6 with 
respect to the y axis from normal incidence (Fig. 3). The 
angles range from 6^ = 0° to6> = 65° with a step AO = 
2.5°. The magnitude of the wavevector k\\ on the sample 
plane is given by k\\ — {uo/c) sin6>, where c is the speed of 
light. 

We perform transmission experiments for two polar- 
izations along the x' axis (parallel configuration) and y 
axis (perpendicular configuration). We denote the elec- 
tric field of the incident wave as E^ and the projection 
of E to the sample plane as E\\. For parallel or perpen- 
dicular configurations, E\\ is parallel or perpendicular, 
respectively, to k\\. Wire-grid polarizers near the emit- 
ter and detector are adjusted so that the emitted and 
detected fields have the same polarization. 
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FIG. 3. (Color online) Schematic view of the transmission 
experiment. The plane of a sample is represented by the co- 
ordinate system (x,y) shown in Fig. 1(a). 
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FIG. 4. (Color online) Parallel configuration (E\\ \\ k\\). 
(a) Experimentally obtained transmission diagram of KBDRs. 
Transmission minima between 0.21 THz and 0.28 THz are ob- 
served and theoretically fitted by a dotted line. The solid 
line corresponds to ^ = 20°. (b) Transmission spectrum for 
= 20° obtained by simulation and experiment, (c) Sur- 
face electric charge distribution obtained by simulation at the 
transmission minimum 0.255 THz for 6 — 20° in the unit cell. 



IV. RESULTS 

Figure 4(a) displays the transmission spectrum for par- 
allel configuration. The wavevectors are estimated as 
k\\ = (oj/c) sin 0. Transmission spectrum minima are ob- 
served from 0.21 THz to 0.28 THz. With an increase of 
wavenumber, the frequency of the transmission minimum 
decreases from 0.28 THz at the F point and approaches 
0.21 THz at the M point. 

For further investigation, we calculate the electromag- 
netic response of KBDRs for parallel configuration. A 
commercial finite-element method solver (Ansoft HFSS) 
is used. In the simulation, a plane THz wave is injected 
into perfectly conducting KBDRs at an incident angle 6. 
By using periodic boundaries with some phase shifts, the 
transmission and the electromagnetic fields in the unit 
cell are calculated for an oblique incident plane wave. 
The measured transmission spectra for = 20° are com- 
pared with the simulation in Fig. 4(b). The frequency 
of the transmission minimum and shape of the curve for 
the simulation are consistent with the experimental re- 
sult, which confirms the validity of the assumption of 
perfect conductors. 

Figure 4(c) shows the calculated distribution of surface 
electric charges at a minimum (0.255 THz) for = 20°. 



This mode corresponds to the middle band. Disks along 
the X axis are alternately charged. The in-phase currents 
along ai and a2 are excited by the electromotive force 
due to E. No resonance appears for ^ = because the 
current flowing into a disk is balanced by the current 
flowing out of it. For the excitation of this mode, a phase- 
shifted field in the x-direction is needed. 

By using Eq. (4), the fitting parameters are obtained 
from experimental data as uJo/{27r) = 0.105 THz and k = 
0.103. The resultant dispersion curve is represented as a 
dotted curve in Fig. 4(a). Positive charges on one disk 
induce negative charges on the other; therefore, a>: < is 
ordinarily expected in the static limit (cj ^ 0). It seems 
strange that Hi would be positive. In our situation, it can 
be explained by a retardation effect. ^^ The phase shift 
between nearest disks is given by (ujc/c) x 1/V3 ~ 0.77 x tt 
at the center frequency co'c/(27r) = 0.25 THz of the middle 
band. The near tt shift leads to a>: > 0. Although hz 
depends on frequency, we can approximately regard it as 
a constant between 0.2 THz and 0.3 THz. 

Figure 5(a) displays the transmission spectrum for per- 
pendicular configuration. Unlike in the case of parallel 
configuration, the flat band of the transmission minima 
is observed around 0.28 THz. 

In order to confirm that the flat band is due to the 
interference of a spoof plasmon, we perform a simula- 
tion for perpendicular configuration. A calculated trans- 
mission spectrum by simulation for = 20° is shown 
in Fig. 5(b) with the experimental data. We see a good 
agreement in the frequency of the transmission minimum 
and the shape of the curve. The calculated distribution 
of surface electric charges at a minimum (0.278 THz) for 
= 20° is shown in Fig. 5(c). The resonant mode has an- 
tisymmetric amplitudes on the right two disks and there 
is no charge stored on the left disk. This mode can be 
constructed by the localized modes shown in Fig. 2(b). 
Therefore, the flat transmission minima are caused by 
the topological nature of the kagome lattice. The mode 
is excited by anti-phase electromotive force caused by E 
along bars parallel to ai and a2 . The electromotive force 
along vertical bars does not contribute to the storage of 
charges on disks because the currents flowing into and 
out of a disk are balanced. In the case of ^ = 0, the 
current flowing into a disk is equal to the current flowing 
out of it and the flat-band mode cannot be excited. 

A dotted line in Fig. 5(a) represents the highest band 
given by Eq. (4) with the previously derived parameters 
cJo/(27r) = 0.105 THz and k = 0.103. It fits well with the 
minima experimentally obtained. The bend of the fiat 
band caused by coupling to second (or higher) nearest 
sites is negligibly small, so the assumption of only nearest 
mutual disk coupling is appropriate. 



V. DISCUSSION 

Our bar-disk resonators (BDRs) are obtained by in- 
verting the metallic area and empty space of the slit-hole 
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FIG. 5. (Color online) Perpendicular configuration [E ± 
fcy). (a) Experimentally obtained transmission diagram 
of KBDRs. A flat transmission band is observed around 
0.28 THz and theoretically fitted by the dotted line. The solid 
line corresponds to ^ = 20°. (b) Transmission spectrum for 
Q — 20° obtained by simulation and experiment, (c) Sur- 
face electric charge distribution obtained by simulation at a 
transmission minimum of 0.278 THz for Q — 20° in the unit 
cell. 



resonators (SHRs)^^'^^ composed of slits and holes en- 
graved on an ultra thin metallic plate. The BDRs and 
SHRs are complementary structures related through the 
Babinet's principle, ^^^^ based on the electric/magnetic 
reciprocity of a vacuum. We denote a pair of an elec- 
tric field El and a magnetic field H as (E^H). Due 
to the Babinet principle, the transmittance of a com- 
plementary metallic screen with an complementary in- 
cident wave {E',H') = {ZqH,-E/Zq) is equal to the 
reflectance of the original metallic screen illuminated by 
an incident wave {E^H)^ where Zq is the impedance of 
a vacuum. Thus, the transmission peaks in SHRs cor- 
respond to the transmission minima in BDRs. This fact 
shows the duality of the Lagrangians of SHRs and BDRs. 

Electromagnetic flat bands for all crystal directions 
have been reported for photonic crystals with square 
symmetry, theoretically^^ and experimentally.^^ In this 
case, the flat bands are formed due to good lateral con- 
finement (high Q factor) of the quadropole modes, which 
lack preferential coupling directions, at defects of pho- 
tonic crystals. On the other hand, the flat band for 
KBDRs is not caused by highly confined modes, but 
by the topological nature of the kagome lattice. The 



kagome lattice prevents spoof plasmons from propagat- 
ing despite the existence of strong coupling. Thus, the 
physical origin of the flat band on KBDRs differs from 
that for photonic crystals. ^^'^^ The flat band for propa- 
gating modes has been theoretically predicted for square 
waveguide networks. ^^ Our system is considered as an 
experimental realization of the flat band for the propa- 
gating mode. 

The flat band in the kagome lattice comes from local 
interference effects. The global symmetry (i.e., period- 
icity of the lattice) is not necessarily required because 
local symmetries can support the localized mode. The 
resonance independent of the incident angle could be ex- 
pected for the metallic structure having localized modes 
with the same resonant frequency, even without period- 
icity. 



VI. CONCLUSION 

In conclusion, we studied theoretically and experimen- 
tally the electromagnetic flat band on a metallic kagome 
lattice. Kagome-type bar-disk resonators were proposed 
to realize the flat band. A dispersion relation composed 
of three bands was theoretically predicted for KBDRs. 
The highest band was flat for all wave vectors. Two bands 
formed by transmission minima depending on the po- 
larization of the incident terahertz beams were observed 
experimentally. One of the bands corresponded to the 
flat band. Theoretical fitting showed good agreement 
for these modes. By simulation, we revealed that the flat 
band was caused by the topological nature of the kagome 
lattice. 

The fiat band can be applicable to slow light, which 
is useful for the control of group velocity, ^"^'^^ high sensi- 
tive sensing, and other applications. In the fiat band, the 
effective mass of photons becomes very heavy and their 
correlation has an important role. Multiphoton correla- 
tion effects in kagome lattices are important in terms of 
fundamental physics and should be studied in the future. 



ACKNOWLEDGMENTS 

The authors would like to thank A. Yao and T. Hik- 
ihara for technical assistance, and S. Endo for fruitful 
discussions. This research was supported in part by 
Grants-in-Aid for Scientific Research No. 22109004 and 
No.22560041, the Global COE program "Photonics and 
Electronics Science and Engineering" at Kyoto Univer- 
sity, the Program for Improvement of Research Environ- 
ment for Young Researchers from the Special Coordina- 
tion Funds (SCF) for Promoting Science and Technol- 
ogy commissioned by the Ministry of Education, Cul- 
ture, Sports, Science and Technology (MEXT) of Japan 
(T.O.), the Research Foundation for Opto- Science and 
Technology (T.O.), and research grants from the Murata 
Science Foundation (T.O. and T.N.). 



* nakata@giga.kuee.kyoto-u.ac.jp 

^ kitano@kuee.kyoto-u.ac.jp 

^ J. L. Atwood, Nature Mater. 1, 91 (2002). 

^ E. H. Lieb, Phys. Rev. Lett. 62, 1201 (1989). 

2 A. Mielke, J. Phys. A 24, L73 (1991). 

^ A. Mielke, J. Phys. A 24, 3311 (1991). 

^ H. Tasaki, Phys. Rev. Lett. 69, 1608 (1992). 

^ S. D. Huber and E. Altman, Phys. Rev. B 82, 184502 

(2010). 
^ G. Moller and N. R. Cooper, Phys. Rev. Lett. 108, 045306 

(2012). 
^ C. Wu, D. Bergman, L. Balents, and S. Das Sarma, Phys. 

Rev. Lett. 99, 070401 (2007). 
^ H. Takeda, T. Takashima, and K. Yoshino, J. Phys. Con- 
dens. Matter 16, 6317 (2004). 
^° S. Endo, T. Oka, and H. Aoki, Phys. Rev. B 81, 113104 

(2010). 
^^ T. Baba, Nature Photon. 2, 465 (2008). 
^^ J. B. Pendry, L. Martin- Moreno, and F. J. Garcia- Vidal, 

Science 305, 847 (2004). 
^^ F. J. Garcia- Vidal, L. Martin- Moreno, and J. B. Pendry, 

J. Opt. A: Pure Appl. Opt. 7, S97 (2005). 
^^ C. R. Williams, S. R. Andrews, S. A. Maier, A. L 
Fernandez-Dominguez, L. Martin- Moreno, and F. J. 
Garcia- Vidal, Nature Photon. 2, 175 (2008). 
^^ S. A. Maier, S. R. Andrews, L. Martin- Moreno, and F. J. 
Garcia- Vidal, Physical Review Letters 97, 176805 (2006). 
^^ V. G. Veselago, Sov. Phys. Usp. 10, 509 (1968). 
^^ R. A. Shelby, D. R. Smith, and S. Schultz, Science 292, 

77 (2001). 
^^ J. Pendry, A. Holden, D. Robbins, and W. Stewart, IEEE 

Trans. Microwave Theory Thech. 47, 2075 (1999). 
^^ J. B. Pendry, Phys. Rev. Lett. 85, 3966 (2000). 
2° U. Leonhardt, Science 312, 1777 (2006). 



^^ J. B. Pendry, D. Schurig, and D. R. Smith, Science 312, 

1780 (2006). 
^^ D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. 

Pendry, A. F. Starr, and D. R. Smith, Science 314, 977 

(2006). 
22 H. Liu, Y. M. Liu, T. Li, S. M. Wang, S. N. Zhu, and 

X. Zhang, Phys. Status Solidi B 246, 1397 (2009). 
2^ N. Biggs, Algebraic Graph Theory, 2nd ed. (Cambridge 

University Press, New York, 1994). 
2^ L. Solymar and E. Shamonina, Waves in Metamaterials 

(Oxford University Press, Oxford, 2009). 
26 H. Liu, T. Li, Q. J. Wang, Z. H. Zhu, S. M. Wang, J. Q. 

Li, S. N. Zhu, Y. Y. Zhu, and X. Zhang, Phys. Rev. B 79, 

024304 (2009). 
2^ C. Zhu, H. Liu, S. M. Wang, T. Li, J. X. Cao, Y. J. Zheng, 

L. Li, Y. Wang, S. N. Zhu, and X. Zhang, Optics Express 

18, 26268 (2010). 
2^ J. Jackson, Classical electrodynamics, 3rd ed. (Wiley, New 

York, 1999). 
^^ M. Born and E. Wolf, Principles of optics, 7th ed. (Perga- 

mon, Oxford, 1999). 
2° F. Falcone, T. Lopetegi, M. A. G. Laso, J. D. Baena, 

J. Bonache, M. Beruete, R. Marques, F. Martin, and 

M. Sorolla, Phys. Rev. Lett. 93, 197401 (2004). 
^^ H. Altug and J. Vuckovic, Appl. Phys. Lett. 84, 161 (2004). 
^2 H. Altug and J. Vuckovic, Appl. Phys. Lett. 86, 111102 

(2005). 
^^ E. Feigenbaum and H. A. Atwater, Phys. Rev. Lett. 104, 

147402 (2010). 

34 Y Tamayama, T. Nakanishi, Y. Wakasa, T. Kanazawa, 
K. Sugiyama, and M. Kitano, Phys. Rev. B 82, 165130 
(2010). 

35 Y^ Tamayama, T. Nakanishi, and M. Kitano, Phys. Rev. 
B 85, 073102 (2012). 



